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Abstract 

We analyze the nonlinear elliptic problem An = rj+uv 3 on a bounded domain ft of M. N with Dirich- 
let boundary conditions. This equation models a simple electrostatic Micro-Electromechanical System 
(MEMS) device consisting of a thin dielectric elastic membrane with boundary supported at above a 
rigid ground plate located at —1. When a voltage -represented here by A- is applied, the membrane 
deflects towards the ground plate and a snap-through may occur when it exceeds a certain critical value 
A* (pull-in voltage). This creates a so-called "pull-in instability" which greatly affects the design of many 
devices. The mathematical model lends to a nonlinear parabolic problem for the dynamic deflection of 
the elastic membrane which will be considered in forthcoming papers II 1 1 and |12|. For now, we focus 
on the stationary equation where the challenge is to estimate A* in terms of material properties of the 
membrane, which can be fabricated with a spatially varying dielectric permittivity profile /. Applying 
analytical and numerical techniques, the existence of A* is established together with rigorous bounds. 
We show the existence of at least one steady-state when A < A* (and when A = A* in dimension TV < 8) 
while none is possible for A > A* . More refined properties of steady states -such as regularity, stability, 
uniqueness, multiplicity, energy estimates and comparison results- are shown to depend on the dimension 
of the ambient space and on the permittivity profile. 



Key words: MEMS; pull-in voltage; power law permittivity profile; minimal solutions. 



Contents 

11 Introduction! 

12 Pull-in Voltage! 



2.1 Lower bounds for A : 



2 i 2_^rjrjeH3mmds i Jbr_^ 



2.3 Numerical estimates for A*l 9 



|3 Min imal Positive Solutio ns! 10 

3.1 Spectral properties of minimal solutionsl 11 



3.3 Existence of solutions at A = A' 



J^^^Enra^-^^s^m^s^fflfljSaUlaj^S 13 



17 

|4 Uniqueness and Multiplicity of Solutionsl 18 
15 Steady-State: Case of Power-Law Profile! 20 



'Partially supported by the Natural Science and Engineering Research Council of Canada. 

tPartially supported by the Natural Science Foundation of P. R. China (10171036) and by a U.B.C. Graduate Fellowship. 



1 Introduction 



Micro-Electromechanical Systems (MEMS) are often used to combine electronics with micro-size mechanical 
devices in the design of various types of microscopic machinery. MEMS devices have therefore become key 
components of many commercial systems, including accelerometers for airbag deployment in automobiles, 
ink jet printer heads, optical switches and chemical sensors and so on (see for example [IS]). The simplicity 
and importance of this technique have inspired numerous researchers to study mathematical models of 
electrostatic-elastic interactions. The mathematical analysis of these systems started in the late 1960s with 
the pioneering work of H. C. Nathanson and his coworkers \TS] who constructed and analyzed a mass-spring 
model of electrostatic actuation, and offered the first theoretical explanation of pull-in instability. At roughly 
the same time, G. I. Taylor [21] studied the electrostatic deflection of two oppositely charged soap films, and 
he predicted that when the applied voltage was increased beyond a certain critical voltage, the two soap films 
would touch together. Since Nathanson and Taylor's seminal work, numerous investigators have analyzed 
and developed mathematical models of electrostatic actuation in attempts to understand further and control 
pull-in instability. An overview of the physical phenomena of the mathematical models associated with the 
rapidly developing field of MEMS technology is given in [201 • 



The key component of many modern MEMS is the simple idealized electrostatic device shown in Fig. ^ 
The upper part of this device consists of a thin and deformable elastic membrane that is held fixed along its 
boundary and which lies above a rigid grounded plate. This elastic membrane is modeled as a dielectric with 
a small but finite thickness. The upper surface of the membrane is coated with a negligibly thin metallic 
conducting film. When a voltage V is applied to the conducting film, the thin dielectric membrane deflects 
towards the bottom plate, and when V is increased beyond a certain critical value V* -known as pull-in 
voltage- the steady-state of the elastic membrane is lost, and proceeds to touchdown or snap through at a 
finite time creating the so-called pull-in instability. 

A mathematical model of the physical phenomena, leading to a partial differential equation for the 
dimcnsionless dynamic deflection of the membrane, was derived and analyzed in 1 and 14'. . In the 
damping-dominated limit, and using a narrow-gap asymptotic analysis, the dimensionless dynamic deflection 
u = u(x, t) of the membrane on a bounded domain O in K 2 , is found to satisfy the following parabolic problem 
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Figure 1: The simple electrostatic MEMS device. 
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for x € O , 



(1.1a) 



u(x, t) = 
u(sg,0) = 



for x E dil , 
for 16(1. 



(1.1b) 
(1.1c) 



An outline of the derivation of (|1.1|) was given in Appendix A of |14j . This initial condition in (|l.lc|) assumes 
that the membrane is initially undeflected and the voltage is suddenly applied to the upper surface of the 
membrane at time t = 0. The parameter A > in (|l.la|l characterizes the relative strength of the electrostatic 



and mechanical forces in the system, and is given in terms of the applied voltage V by 



x -^r^r> (L2) 

where d is the undeflected gap size (see Fig.^), L is the length scale of the membrane, T e is the tension of the 
membrane, and £o is the permittivity of free space in the gap between the membrane and the bottom plate. 
In view of relation (jl.2Jl , we shall use from now on the parameter A and A* to represent the applied voltage 
V and pull-in voltage V* , respectively. Referred to as the permittivity profile, f(x) in (1.1a) is defined by 
the ratio 

f(x) = ^, (1.3) 

where £2(2;) is the dielectric permittivity of the thin membrane. 

There are several issues that must be considered in the actual design of MEMS devices. Typically one 
of the primary goals is to achieve the maximum possible stable deflection before touchdown occurs, which is 
referred to as pull-in distance (cf. ^1] and ^H]). Another consideration is to increase the stable operating 
range of the device by improving the pull-in voltage A* subject to the constraint that the range of the 
applied voltage is limited by the available power supply. Such improvements in the stable operating range 
is important for the design of certain MEMS devices such as microresonators. One way of achieving larger 
values of A*, while simultaneously increasing the pull-in distance, is to use a voltage control scheme imposed 
by an external circuit in which the device is placed (cf. This approach leads to a nonlocal problem 

for the dynamic deflection of the membrane. A different approach studied in ^2] and ^2] ^ s to introduce a 
spatially varying dielectric permittivity £2(2;) of the membrane. The idea is to locate the region where the 
membrane deflection would normally be largest under a spatially uniform permittivity, and then make sure 
that a new dielectric permittivity £2(2;) is largest -and consequently the profile f(x) smallest- in that region. 

This latter approach requires the membrane having varying dielectric properties, a framework investigated 
recently in JHj and 52- ^ n El J- Pelesko studied the steady-states of (1.1), when f(x) is assumed to be 
bounded away from zero, i.e., 

< C*o < f(x) < 1 for all x e O. (1.4) 

He established in this case an upper bound Ai for A*, and derived numerical results for the power-law 
permittivity profile, from which the larger pull-in voltage and thereby the larger pull-in distance, the existence 
and multiplicity of the steady-states were observed. Recently, Y. Guo, Z. Pan and M. Ward studied in 
the dynamic behavior of (1.1), which is also of great practical interest. They considered a more general class 
of profiles j{x), where the membrane is allowed to be perfectly conducting, i.e., 

0</(x)<l for all x£Q, (1.5) 

with f(x) > on a subset of positive measure of £1. By using both analytical and numerical techniques, 
they obtained larger pull-in voltage A* and larger pull-in distance for different classes of varying permittivity 
profiles. They obtained new upper bounds A2 for the pull-in voltage A* that correspond to the more general 
profiles f(x) satisfying (|1.5|l . Moreover, having estimated A* numerically as some saddle-node bifurcation 
value, they showed that A* is generally strictly smaller than both Ai and A2 (cf. Table 1 of (11)1. 

In this paper, we shall focus on the stationary deflection of the elastic membrane, leaving the dynamic 
case to our forthcoming papers JO and |12| . For convenience, we shall set v = —u in such a way that our 
discussion will center on the following elliptic problem 

< v < 1 x € 0; 1 M 

v = x e dQ. 

We shall continue the investigation of optimal upper and lower bounds for the pull-in voltage, and how 
they relate to the permittivity profile / which will be assumed to satisfy l|1.5J) throughout the paper unless 
mentioned otherwise. We shall also discuss the issues of existence, multiplicity, and other related properties 
of steady-states for (S)\, and their remarkable dependence on space-dimension. 



This paper is organized as follows: In §2 we mainly show the existence of a specific pull-in voltage and 
establish lower and upper bound estimates. For that, we shall write cu N for the volume of the unit ball fli(O) 
in M. N , and for any bounded domain T in M. n , we denote by /x r the first eigenvalue of —A on ^(I^and by 
4> r (resp., tp r ) the corresponding positive eigenfunction normalized with J T (f> r dx = 1 (resp., sup^gp ip r = 1). 
We shall also associate to any domain f2 in R w the following parameter: 

v n = sup{M r -ff(inf ipr); r domain of R N , T D (1.6) 
where H is the function Hit) — t ( t+1 + 2 v^) _ 

v ' (t+l + Vt) 

We then prove the following lower estimates, the upper ones having been established in ^] and |14j . 
Theorem 1.1. There exists a finite pull-in voltage A* > such that 

1- If < A < A*, there exists at least one solution for (S)\; 

2. If X > X* , there is no solution for (S)\. 
Moreover, we have the bounds 

maXi ^rp7R'27sup/(x)VM^ 1 - A - mm \27 inf f(x)'3Lf<l> n dxJ 

Furthermore, if f{x) = \x\ a on Q with a > 0, then we have the more refined lower bound 

X c (a) := [—) < A (1.7) 

In §2.3 we give some numerical estimates on A* to compare them with analytic bounds given in Theorem 
1.1. Note that the upper bound Ai = ( inf^go f( x YJ is relevant only when / is bounded away from 

0, while the upper bound A2 = ^f 2 -^ fnf&n dxj is valid for all permittivity profiles. In the case of a 
uniform permittivity profile / = 1 on Q, where f2 is a strictly star-shaped domain containing 0, we give a 
more explicit upper bound A3 of A* in Proposition 2.4. In particular, we show that A3 = is an upper 

bound in the case where the domain is the unit ball £1 = -Bi(O) C K . 

The issues of uniqueness and multiplicity of solutions for (S)\ with < A < A*, and even mere existence 
for (S)\* seem to be quite interesting. We address these problems beginning in section §3 by first considering 
minimal (positive) solutions of (S)\ defined as follows. 

Definition 1.1. A solution < u\(x) < 1 is said to be a minimal (positive) solution of {S)\, if for any 
solution < u{x) < 1 of (S)\ we have u\(x) < u{x) in fl. 

Our main results in this direction can be stated as follows. 

Theorem 1.2. Under the above assumptions, and with X* as defined in Theorem 1.1, there exists for any 
X < X* , a unique minimal positive classical solution U\(x) of (S)\. It is obtained as the limit of the sequence 
{it„(A; x)} constructed recursively as follows: u a = in Q and for each n > 1, 

(1 ~ Un-i) 2 

Q< u n < 1, x en; C 1 - 8 ) 

u n = , x e dfl. 

Moreover, minimal solutions satisfy the following properties: 

1. For each x £ O, the function X — > u\(x) is strictly increasing and differ entiable on (0, A*); 



2. If 1 < N < 8, then there exists a constant < C(N) < 1 such that || u 



a llc(n) 



< C(N) for all\< A*. 



We refer to Lemma 3.6 in §3.2 for a more general version of Theorem 1.2(2). Based on the results of 
Theorem 1.2, the existence and related properties of minimal solutions at critical voltage A = A* will be 
studied in §3.3. More precisely, we shall establish the following. 

Theorem 1.3. If 1 < N < 8 then u x „ = lim AXA » u x exists in the topology ofC 2 ' a {il) with < a < 1, and 
u x , is the unique classical solution of (S)a* • 

§4 is devoted to the uniqueness and multiplicity of solutions which remarkably depend again on the 
space-dimension. 

Theorem 1.4. Under the above assumptions, with A* defined as in Theorem 1.1, we have: 

1. If N > 2, then for any M > there exists a voltage < \\{M) < A* such that for every A G 
(0, \\{M)), there exists a unique positive solution for {S)\ -namely the minimal solution u\- that 



2. If 1 < N < 8 then there exists < < A* such that (S)\ has at least two solutions for A € (A2, A*). 

A uniqueness result in the spirit of (1) also holds for dimension 1 (resp., dimension 2) with N/2 replaced 
by 1 (resp., 1 + e). However, in spite of above results, issues of uniqueness, multiplicity and other qualitative 
properties of the solutions for (S)\ are still far from being well understood. For example, we conjecture 
that no solution exists for (S*)a* with N > 8 -at least when / = 1. In §5 we shall present some numerical 
evidences for various conjectures relating to the case of power-law permittivity profile f(x) = \x\ a defined in 

a unit ball. It looks like there are two critical exponents a* = -\ + ^27/2 and a**(N) = 4 ~ 6jv+3 ^ (jv ~ 2) 
(which is relevant for N > 8) such that the following four regimes are possible: 

1. There exist exactly two solutions for < A < A*, and one solution for A = A*. This regime occurs 
when N = 1 and a < 1. 

2. There exists exactly one solution for < A < A*, exactly 2 solutions for A^ < A < A* and exactly one 
at A = A*. This regime occurs when N = 1 and 1 < a < a*. 

3. There exists exactly one solution for < A < A^, exactly two solutions for A2 < A < A*, while multiple 
solutions can be obtained for Aj < A < A2. Moreover, the multiplicity becomes arbitrarily large as A 
approaches another critical value A* € (A*, A2), at which there is a touchdown (quenching) solution u 
characterized with || u ||oo= 1- This regime occurs when 

• 2 < N < 7 and a > 0; 

• N > 8 and a** < a. 

4. There is exactly one solution if < A < A* and none for A > A*. This regime occurs when N > 8, and 



We finally mention that the above results can be extended to more general elliptic problems of the form 



satisfies | 




\~alx < M; 



0<a<a 



-Av 



AM-'' 

{i-vY 



v(x) = 0, 



xGdn 



with (3 > 0. Here the critical dimension depends on the parameter (3, and this is the subject of a work in 
progress. 



2 Pull-in Voltage 

In this section, we study the steady-state deflection u which satisfies (S)\, and we establish the existence 
and some estimates on the pull-in voltage A* for (S)\ defined as: 

A* = sup{A > | (S)\ possesses at least one solution} . (2.1) 

In other words, A* is called pull-in voltage if there exist uncollapsed states for < A < A* while there are 
none of them for A > A*. 

Theorem 2.1. There exists a finite pull-in voltage A* > such that 

1. If X < X* , there exists at least one solution for (S)\; 

2. If X > X* , there is no solution for (S)\. 

Moreover, with v n defined by (1.6), we have the lower bound 

^(sup/Or))" 1 < A*. (2.2) 

Proof: We need to show that (S)\ has at least one solution when A < ^ si (sup f(x))~ 1 . Indeed, it is clear 
that u = is a sub-solution of (S) \ for all A > 0. To construct a super-solution of (S)\, we consider a 
bounded domain T D H with smooth boundary, and let (fi r ,ip r ) be its first eigenpair normalized in such a 
way that 

supVr(a;) = 1 and inf ipr(%) ■= s\ > 0. 
xer x ^ n 

We construct a super-solution in the form ?/> = Aip r where A is a scalar to be chosen later. First, we must 
have Aip T > on <9f2 and < 1 — Atp r < 1 in fi, which requires that 

0<a<l. (2.3) 

We also require 



which can be satisfied as long as: 



or 

Asup/(x) <f3(A,T) :=/i r inf{ ff (*A);se[si(r),l]}, (2.6) 

Q 

where g(s) = s(l — s) 2 . In other words, A*sup/(x) > sup{(3(A, T): < a < 1,T D 0}, and therefore it 

o 

remains to show that 

u a = sup{f3(A,T); 0<o<l,rDfl}. (2.7) 

For that, we note first that 

inf g(As) =ram{g(A Sl ),g(A)}. 

s£[si,l] 

We also have that g(As\) < g(A) if and only if A 2 (s\ — 1) — 2A(s\ — 1) + (si — 1) < which happens if and 
only if A 2 (s\ + si + 1) - 2A(s 1 + 1) + 1 > or if and only if cither A < A- or A > A + where 

_ S 1 + 1 + ygT _ 1 _ Si + 1 - ^~S~l _ 1 

+ ~~ s\ + 1 + si ~~ si + 1 - y/sl ' ~ ~ sj + 1 + si ~ si + 1 + y/sl 
Since A_ < 1 < a + , we get that 

f g(Asi) if < A < A_ , 
G(A) = inf , 9 (As) = ^ (2.8) 
^I*!' 1 ! ( 5 (A) if A_ < A< 1. 



We now have that = g'(As x )si > for all < A < A-. And since A_ > ±, we have = g'(A) < for 
all A- < A < 1. It follows that 

sup inf g(As) = sup G(A) = G(A-) = g{AJ) 

0<a<l se[si,l] 0<a<l 

= I (1 I ) 2 

= Si(si + l + 2^T) 

(si + i + V^T) 3 
= H(M n #)) 

which proves our lower estimate. 

Now that we know that A* > 0, pick A € (0, A*) and use the definition of A* to find a A G (A, A*) such 
that (S)x has a solution u x , i.e, 



A f(x ) 

-Au- X = Jy > x<Efl; u- x = 0, xedn. 



and in particular — Au x > ^I^y for x € fl which then implies that u x is a super-solution of (S)\. Since 
u = is a sub-solution of (S)\, then we can again conclude that there is a solution u\ of [S)\ for every 
A € (0,A*). 

It is also easy to show that A* is finite, since if (S)\ has at least one solution < u < 1, then by 
integrating against the first (positive) eigenfunction <p n , we get 

+00 > fj, n > fi a [ u<p a = - [ uA(f> n = - [ <f) Q Au = A [ ,^"\ 2 dx >\[ <j) n fdx (2.9) 
Jn Jn Jn Jn K 1 ~ u ) Jn 

and therefore A* < +oo. The definition of A* implies that there is no solution of (S)\ for any A > A*. ■ 
2.1 Lower bounds for A* 

It is desirable to seek more computationally accessible lower bounds on A*. For that we consider for every 
subset T C M. N and any function / on T such that < / < 1, the corresponding pull-in voltage A*(r, /), 
that is the value A* defined above for the problem 

- Au = - y ' xeT, 2.10a 
(1 — u) z 

< u < 1 x £ r, (2.10b) 

u = xedT. (2.10c) 

We need the following result which can be found in [2] (Theorem 4.1). 
Lemma 2.2. For any bounded domain T in 18L N and any function f on T such that < / < 1, we have 

X*(TJ)>X*(B R ,r) 

where Br = Bji(0) is the Euclidean ball in R N with radius R > and with volume \Br\ = |r|, and where 
f* is the Schwarz symmetrization of f . 

We now establish the following refined lower bounds for A* of (S)\. 

Lemma 2.3. We have the following lower bound for A*. 

Moreover, if f(x) — \x\ a on fi with a>0, then we have 

A ^4(2 + q)(iV + a)^^ (2.12) 



27 \\Q 



Proof: Setting R = ( — J , it suffices -in view of Lemma 2.2- and since sup /* = sup / , to show that 

\ U N / Rn o 



8JV 

A* > ^- , (2.13) 

-27i? 2 su Po /' 1 ; 



for the case where fl = Br. In fact, the function w{x) = |(1 — ^r) satisfies on Br 

2N 2N{\-\) 2 1 



-Aw 



3R 2 3R 2 (I-3) 2 

> 87V f(x) 



27i? 2 su Po /[i_ 
87V f{x) 



21R 2 su Vn f {\-w) 2 • 

So for A < 27 sap — J i w 1S & super-solution of (S)\ in Br. Since on the other hand wq = is a sub-solution 
of (S)x and wq < w in Br, then there exists a solution of (S)\ in Br which proves l|2.13[l and hence H2.11fl . 

In order to prove (|2.12|l . it suffices to note that w(x) = A(l — R2+ a ) is a super- solution for (S)\ on Br 
provided A < 4 ^ 2 ^]a+l +a ' ) ■ This completes the proof of Lemma 2.3. ■ 



2.2 Upper bounds for A* 

We note that (|2.9|l already yields a finite upper bound for A*. However, Pohozaev-type arguments such as 
the one used in |14j can be used to establish better and more computable upper bounds. In this subsection, 
we establish these estimates and hence complete the proof of Theorem 1.1. 

We shall consider problem (S)\ in the case where f2 C is a strictly star-shaped domain containing 0, 
meaning that f2 satisfies the additional property that there exists a positive constant a such that 

x ■ v > a / ds for all x S d£l , (2-14) 
Jon 

where v is the unit outer normal to £1 at x £ dQ. 

Proposition 2.4. Suppose f = 1 and that the strictly star-shaped domain f2 C M. N satisfies \2.14\ - Then 
the pull-in voltage A*(0) satisfies: 

(N + 2) 2 

A *^ A -Wr (2 - 15) 

In particular, if ft = B\ (0) C R w then we have the bound 

(N + 2) 2 
A*(S!(0))<^^. 

Proof: Recall the well-known Pohozaev's identity: If u is a solution of 

Au + Xg(u) = for ie!l, 
u = for x 6 dil , 

then 

NX f G(u)dx- A / ug{u)dx=\ [ ( x ■ v )(^-f ds , (2.16) 



where G(u) — J? g(s)ds. Applying it with g(u) = 7j—^p and G{u) — jzt^, it yields 



A f u(N + 2 - 2Nu) _ 1 
(1 - uf ~ 2 



an 



> 



dv' 

2 



(2.17) 



= -( - / Audx) 



aA^ 
2 



dx 



{l-uf 

where we have used the Divergence Theorem and Holder's inequality 



Since 



du 

'an dv 
u(N + 2 - 2Nu) 



N + 2,2 {N + 2) 2 



1/2 



an 



ds 



an 



1/2 



(1-U) 2 



dx = 



< 



- 2N(u- 
{N + 2) 2 



8N 



AN 

d/x 



8N 



(i- u y 



-dx 



we deduce from (|2.17|l that 



(N + 2f 
8N 

which implies the upper bound (|2.15|) for A*. 



> aA 



(1 - u) 2 
dx 



{1-uf 



> a\\n\ , 



Finally, for the special case where SI = -Bi(O) C M. N , we have a = — with ui N = |Bi(0)|, and hence the 



bound A* < A 3 = 



For a general domain SI, the following upper bounds on A* (SI) established in and ^2] respectively, 
complete the proof of Theorem 1.1. 



Proposition 2.5. (1) If f satisfies < Co < f(x) < 1 on SI, then 

X " Al= 27Cb- 

(2) If f satisfies < f{x) < 1 on SI, and if f > on a sei of positive measure, then 

-l 



A, < A 2 = ^ 



f(j> a dx 



(2.18) 



(2.19) 



-ffere fi n and <p n are the first eigenpair of —A on H^(Vl) with J^^^dx = 1. 

2.3 Numerical estimates for A* 

In the computations below we shall consider two choices for the domain SI, 

SI : [—1/2, 1/2] (Slab); SI : x 2 + y 2 < 1 (Unit Disk) . 



(2.20) 



To compute the bounds Ai and A2, we must calculate the first eigenpair ji n and <j) a of — A on SI, normalized 
by J n 4> n dx — 1, for each of these domains. A simple calculation yields that 



= z o ~ 5 - 783 i 



7T \ X + 



(Slab) ; 



Ji(zo) 



Jo( z o\ x \) 7 (Unit Disk) . 



(2.21a) 
(2.21b) 



Exponential Profiles: 



ft 


a 


A 


A* 


Ai 


A 2 


(Slab) 





1.185 


1.401 


1.462 


3.290 


(Slab) 


1.0 


1.185 


1.733 


1.878 


4.023 


(Slab) 


3.0 


1.185 


2.637 


3.095 


5.965 


(Slab) 


6.0 


1.185 


4.848 


6.553 


10.50 


(Unit Disk) 





0.593 


0.789 


0.857 


1.928 


(Unit Disk) 


0.5 


0.593 


1.153 


1.413 


2.706 


(Unit Disk) 


1.0 


0.593 


1.661 


2.329 


3.746 


(Unit Disk) 


3.0 


0.593 


6.091 


17.21 


11.86 



Table 1: Numerical values for pull-in voltage A* with the bounds given in Theorem 1.1. Here the exponential 
permittivity profile is chosen as $2.2ty) . 



Power-Law Profiles: 



n 


a 


A c (a) 


A* 


Ai 


A 2 


(Slab) 





1.185 


1.401 


1.462 


3.290 


(Slab) 


1.0 


3.556 


4.388 


oo 


9.044 


(Slab) 


3.0 


11.851 


15.189 


oo 


28.247 


(Slab) 


6.0 


33.185 


43.087 


oo 


76.608 


(Unit Disk) 





0.593 


0.789 


0.857 


1.928 


(Unit Disk) 


1.0 


1.333 


1.775 


oo 


3.019 


(Unit Disk) 


5.0 


7.259 


9.676 


oo 


15.82 


(Unit Disk) 


20 


71.70 


95.66 


oo 


161.54 



Table 2: Numerical values for pull-in voltage A* with the bounds given in Theorem 1.1. Here the power-law 
permittivity profile is chosen as 

Here Jo and J\ are Bessel functions of the first kind, and z$ « 2.4048 is the first zero of Jq(z). The bounds 
Ai and A 2 can be evaluated by substituting (|2.21() into (|2.18() and i|2.19ll . Notice that A 2 is, in general, 
determined only up to a numerical quadrature. 

Using Newton's method and COLSYS [Q, one can also solve the boundary value problem (S)\ and 
numerically calculate A* as the saddle-node point for the following two choices of the permittivity profile: 

(Slab) : f(x) = \2x\ a , (power-law) ; f{x) = e a ^ 2 - 1/4) (exponential) , (2.22a) 
(Unit Disk) : f{x) = \x\ a , (power-law) ; f(x) = e"^! 2 - 1 ) , (exponential) , (2.22b) 

where a > 0. Table 1 contains numerical values for A* in the case of exponential profiles, while Table 2 deals 
with power-law profiles. What is remarkable is that Ai and A 2 are not comparable even when / is bounded 
away from and that neither one of them provides the optimal value for A*. This leads us to conjecture 
that there should be a better estimate for A*, one involving the distribution of / in ft, as opposed to the 
infimum or its average against the first eigenfunction <j) n . 

3 Minimal Positive Solutions 

In this section, we are concerned with minimal positive solutions for {S)\. We establish their existence, 
uniqueness and other related properties. We consider the case A G (0, A*) in §3.1, and A = A* in §3.3, but 
first we give a recursive scheme for the construction of such solutions. 

Theorem 3.1. For any A < A* there exists a unique minimal positive solution u\ for (S)\. It is obtained 



as the limit of the sequence {u n (X;x)} constructed recursively as follows: uq = in fi and for each n > 1, 

-Au n — A/jk_ a;€fi; 

(1 - U„_l) 2 
Q< u n < 1, x S fi; 

u„ = , x e 9fi. 

Proof: Let u be any positive solution for (S)a, and consider the sequence {u n (X; x)} defined in Clearly 
u(x) > uq = in fi, and whenever u(a;) > u n -i in 0, then 

- A(u - u n ) = Xf(x) [^-L-. - (i _ 1_ ] > o , xefi 
w — Un = , a; € 957 . 

The maximum principle and an immediate induction yield that 1 > u{x) > u n in fi for all n > 0. In a 
similar way, the maximum principle implies that the sequence {u n (X;x)} is monotone increasing. Therefore, 
{it n (A; a;)} converges uniformly to a positive solution ua(x), satisfying u(x) > ttx(a;) in fi, which is a minimal 
positive solution of (S')a- It is also clear that u\(x) is unique. ■ 

Remark 3.1. Let g(x, £, fi) be the Green's function of Laplace operator, with g(x, £, fi) = on dfi. Then the 
iteration in (|3.1|) can be replaced by: uo = in fi and for each n > 1, 

MA^) = Ay n(1 _ Un i(A;0)2d c, *En ; (3<2) 

u„(A; ar) = , a; € <90 . 

The same reasoning as above yields that lirun^oo u n (X;x) = u\(x) for all x G fi. 

The above construction of solutions yields the following monotonicity result for the pull-in voltage. 

Proposition 3.2. If fix C fi2, i/ien A*(fix) > A*(fi2) anrf the corresponding minimal solutions satisfy 
u ni (A, x) < u Q2 (A, x) on fix for every < A < A* (O2) . 

Proof: Again the method of sub/super-solutions immediatly yields that A*(fix) > A*(fi 2 ). Now consider for 
i = 1,2, the sequences {u n (X, x, fi^)} on fij defined by (|3.2[1 where g(a;, £, fij) are the corresponding Green's 
functions on fij. Since fix C fi2, we have that g(a;,£,fix) < <?(a;,£, f^) on fix. Hence, it follows that 

Ul (A, x, n 2 ) = a / /(Ogfo £, n 2 R > a / /(Os(z, C, fiiR = «i(a, x, no 



on fix. By induction we conclude that it„(A, x, fi 2 ) > u n (X, x, fix) on fix for ft- On the other hand, since 
u n (X,x, fi2) < u n +x(A, a;, fi 2 ) on fi 2 for n, we get that u n (A,a;,fix) < u n (A, a;) on fix, and we are done. ■ 



3.1 Spectral properties of minimal solutions 



L u ,x = -A - - J (3.3) 



For a further study of minimal (positive) solutions, we now consider for each positive solution u of (S)\, the 
operator 

2Xf_ 

(1 - uf 

associated to the linearized problem around u. We see that minimal solutions in the above sense correspond 
to variational solutions that are local minimizers. We denote by /i(A, u) the smallest eigenvalue of L u ^\, that 
is the one corresponding to the following Dirichlet eigenvalue problem 

-A4>-^^4, = t i(X,u)cf>, xefi; (3.4a) 

(1 — 7i) J 

= xedfl. (3.4b) 



In other words, 

H{X, u) = mf r —— . 

Proposition 3.3. The following hold: 

1. X* = sup{A; L u t x has positive first eigenvalue for minimal solution u\ of (S)\\. 

2. If < A < A* then the smallest eigenvalue fix := fi(X,u\) of L Ux .\ -corresponding to the minimal 
solution u\- is positive and X — > fx\ is decreasing on (0, A*). 

For Proposition 3.3, we need the following crucial lemma. 

Lemma 3.4. Suppose u is a positive solution of (S)\, and let /i(A, u) be the corresponding first eigenvalue. 
Consider any -classical- supersolution v of (S)\, that is 



Av > n J v "; 2 xen, (3.5a) 



Xf{x) 
(1-v) 

< v{x) < 1 xen (3.5b) 

v = ieffl. (3.5c) 



If /i(A, u) > then v > u on Q, and if fJ,(X, u) = then v = u on Q. 

Proof: For a given A and x G O, use the fact that f(x) > and that t — > ^z^s i s convex on (0, 1), to 
obtain 

xm 

[1 - (u + t(v - u))} 2 

for t G [0, 1]. Note that (|3.6|l is an identity at r = 0, which means that the first derivative of the left side 
for H3.6fl with respect to r is nonnegative at r = 0, i.e, 



-A(u + t{v-u))-- — - ^ > o xen, (3.6) 



2Xf(x) 
{1-u) 3 

v-u = xedfl. (3.7b) 



A(v-u)- n \.-! 3 (v-u) >0 ie!l, (3.7a) 



Thus, the maximal principle implies that if fi(X, u) > we have v > u on f2, while if /i(u) =0we have 



2A/(x) 
(1 - u) 3 



-A(v-u)- £ J \' (v-u)=0 xeft. (3.8) 



In the latter case the second derivative of the left side for H3.6JI with respect to r is nonnegative a t = 
again, i.e, 

6Xf( ^.{v-u) 2 >0 xe n : (3.9) 



(1 — u) 4 

From (|3.9|l we deduce that v = u in £1 \ CIq, where 

Qa = {x G n : f{x) = for x G f2} . (3.10) 

On the other hand, i|3.8[l reduces to 

—A(v - u) = x G £!o , 
i) — u = a; G <9£!o , 

which implies v = u on fio- Hence if /i(A, u) = then v = it on f2, which completes the proof of Lemma 3.4. 
■ 

Proof of Proposition 3.3: (1) Let 

A** = sup{A; L u \ has positive first eigenvalue for minimal solution u\ of (S)\}. 



It is clear that A** < A* , so it suffices to prove that there is no minimal solution for (S)^ with /i > A** . In 
fact, suppose w is a minimal solution of (S)x**+s with 5 > 0, then we would have for A < A**, 

(i — wy (i — u?)^ 

Since the minimal solutions u\ satisfy — Au x = ^l^yi x £ Q for all < A < A**, it follows from Lemma 

3.4 that 1 > w > u x for all < A < A**. Consequently, u = \im\y\** u x would exist. Now from the 
definition of A** and Lemma 3.4, we must have w = u and 5 = on f2 which is a contradiction, and hence 
A** = A*. 

(2) From the first part we conclude that if < A < A* and u = u x , then the smallest eigenvalue of 
—A — ^zffl is positive. Applying the maximum principle, it is easy to show that u x (x) is increasing with 
respect to A (More details can be found in the proof of Theorem 1.2(1) below). That fi x is decreasing with 
respect to A follows now easily from the variational characterization of ^i x and the convexity of (1 — u)~ 3 
with respect to u. ■ 

Remark 3.2. For the case where f(x) > on f2, Lemma 3 of .?. gives ^(1,0) as an upper bound for A** (= A*). 
It is worth noting that our upper bound A in Theorem 1.1 gives a better estimate. Indeed, if / = 1 then 
/i(l, 0) = ^tfi/2 while the estimate in Theorem 1.1 gives -^p- for an upper bound. 

Proof of Theorem 1.2(1): By Theorem 3.1, it suffices to prove that for each x E £1, the function A — * u\(x) 
is differentiable and strictly increasing on (0, A*). Setting F(\,u\(x)) — —Au\ — , Proposition 3.3 

then implies that F Ux (A, u x ) on f2 is invertible for < A < A*. It then follows from the Implicit Function 
Theorem that u x (x) is differentiable with respect to A. 

Consider now for Ai < A2 < A*, their corresponding minimal positive solutions u\ 1 and u\ 2 and let u* be 
a positive solution for (S)\ 2 . For the monotone increasing series {u n (Xi]x)} defined in l|3.1l) . we then have 
u* > uq(Ai; x) = 0, and if w n _i(Ai; x) < u* in f2, then 

- A(u* - u n ) = f(x) [ A2 - Al ] > , xeil 
1 (l — u*y (i — " 



, xedn. 



So we have u„(Ai; x) < u* in fl. Therefore, u\ 1 — linin^oo u n (Xi;x) < u* in f2, and in particular u\ 1 < u\ 2 



in n. Therefore, > for all x e fl. 

Finally, by differentiating (S)\ with respect to A we get 

du x 2A/(x) du x _ f(x) 
dX (1-uJ 3 dX (l-u x f~ ' ^ 

^ > , xedn. 

dX 

Applying the strong maximum principle, we conclude that > on f2 for all < A < A*. 

3.2 Energy estimates and regularity 

We start with the following easy observation. 

Lemma 3.5. Any positive (weak) solution u in H^fl) of (S)x satisfies J n n _ M >2 dx < cxd. 
Proof: Since u € Hq(£1) is a positive solution of (S)x, we have 

which implies that 

/ 



Lj^ iC+ b c ^ + i'^ c+Cc Lw 



uf 



with e > 0. Therefore, by choosing e > small enough, we conclude that J Q ^l u yi < oo . ■ 

That f/(l — u) E L 2 (£l) is unfortunately not sufficient to obtain regularity results for the solutions. 
However, we now show that the situation is much better if //(l — u) has better integrability properties. 

Theorem 3.6. For any bounded domain Q C M. N and any constant C > there exists < K (C, N) < 1 
such that a positive weak solution u of (S)\ (0 < A < A*) is a classical solution and || u || c(n) < K(C,N) 
provided one of the following conditions holds: 

1. N=land\\ J ^ w \\ Ll(n) <C. 

2. N = 2 and WjT^\\ Ll+Hn} < C for some e > 0. 

3. N > 2 and Utt^H w/2 < C. 

Proof: We prove this lemma by considering the following three cases separately: 

(1) If N = 1, then for any I > we write using the Sobolev inequality with constant K(l) > 0, 

K(l) || (1 -„)-i- ly^ < f IvKl-u)- 1 -!]] 2 

Jn 

<CI + C [ f(l-u)~ 5 

J{(l-«)- 3 >/} 

<CI + C [ 8/(1- u)- 2 

J {{l-u)- 3 >I} 

+C f /[(l- M )- 3 + 2(l-,i)- 2 +4(l- U )- 1 ][(l- U )- 1 -l] 2 

J{(l-u)- 3 >/} 

< CI + C + C || (1 - I*)"! - 1 |||oo ({(1 _ u) -3 >/}) / 

^ ' - " J{(l-u)-3>I} U - M J 

<CI + C + Ce(I) || (1 - u)- 1 - 1 ||i_ 

(3.11) 

with s(I) = [ - — ^ . From the assumption f / (1 — u) 3 G ^(Q), we have e(I) — > as 7 — > oo. 
A(i-«)" 3 >-f} ( 2 ~ M ) 

We now choose / such that e{I) < ^£r-, so that the above estimates imply that || (1 — tt) _1 — 1 ||l°°< K(C) . 
Standard regularity theory for elliptic problems now imply that 1/(1 — u) e C 2,a {Vl). Therefore, u is classical 
and there exists a constant K(C, N) which can be taken strictly less than 1 such that || u || < K (C, N) < 1. 

(2) The case when TV = 2 is similar as one can use that Hg embeds in L p for any p < +oo. 

(3) The case when N > 2 is more elaborate and we first show that (1 — u)^ 1 e L q (Q) for all q e (l,oo). 
Since u € -ffg(^) is a solution of (S)\ 1 we already have J < ^- Now we proceed by iteration to show 

that if J n (1 _J )2 +2e < C for some 9 > 0, then J n (1 _ u) l, {1+e) < C. 

Indeed, for any constant 9 > and I > we choose a test function <f> — [(1 — u)~ 3 — 1] min{(l — u)~ 2e , £ 2 }. 
By applying this test function to both sides of (S)\, we have 

A [ f(l-u)- 2 [(l-u)- 3 -l]mm{(l-u)- 2e ,£ 2 }= f Vu • V [((1 - u)' 3 - l) min{(l - u)- 20 , £ 2 }] 
Jn Jn 

= 3 / \Vu\ 2 {l-u)- 4 mm{(l-u)- 2e ,£ 2 } + 29 f |Vu| 2 (l - n)- 20 ' 1 ^ - u)~ 3 - 1] . 

JO J{(1-m)-«<£} 

(3.12) 



ivufa-to-"-"^-!*)- 1 -!] 1 



We now suppose ^i_Jyi+is < C We then obtain from (|3.12l) and the fact that ^J^s < C t jY^u^ij^n ~ ^-Y 
for (1 - u)~ 3 >I>1 that 

|V[((1 - u)- 1 - 1) min{(l - u)- e ,£}}\ 2 

<2 [ \Vu\ 2 (l-u)-' l mm{(l-u)- 20 ,£ 2 } + 29 2 [ 

Jn J{(i-u)- 

= 2 f \Vu\ 2 (l-u)- A mm{(l-u)- 2e ,£ 2 } 
Jn 

+29 2 [ |Vw| 2 (l - u)- 20 - 1 [(1 - u)- 3 - 1 + 1 + (1 - u)- 1 - 2(1 - u)- 2 ] 

J{(i-u)- e <e} 

<CX f f(l-u)- 2 [{l- U y 3 - l]min{(l-u)- 29 ,^ 2 } 
Jn 

<CX [ f(l- u)- 5 mm{(l- u)- 20 ,e 2 } 



(3.13) 



<CI + C 
<CI + C 
<CI + C 



{(!-«)-»>/} 



/(l- M )- 5 min{(l-w)- 2S ,£ 2 } 



/(l - M )~ 3 [(1 - u)- 1 - 1] min{(l - u)- 26 , £ 2 } 



{(l-tt)-3>/} 



((1-li) 3 ) 



([(1 - w)" 1 - l] min{(l - u) -8 ,*}) 

'{(l-u)- s >/> 

<C/ + Ce(I) f |V[((l-u) _1 -l)min{(l-u)-V}]| 



2A r 
JV-2 



with 



e (j) 



(fl -u) 3 ) 



{(l- u )-3>/} v(l - u) 



From the assumption //(l — w) 3 € L™ (f2) we have e(J) — * as / — * oo. We now choose I such that 
e(I) = and the above estimates imply that 



f |V[(1 - u)- 6 - 1 - (1 - u)- e ] | 2 <CI, 



I {(l-u)- e <l} 

where the bound is uniform with respect to £. This estimate leads to 

,2 



(e+i)' 2 f{{i-u)-"<e} I^K 1 ~ M ) 9 X ] 



'{(i-«)- e <£} 
< CJ + C 



26-4 



(1-u) 



-20-3 



< CJ + 



{(i-«)- 9 <^> 

[Ce(l-u)- 29 - 4 + C/e]|Vu| 



{(!-«)-«<«} 



< CI + Ce 

with £ > 0. This means that for e > sufficiently small 

\-0-H 2 _ 



{(l-u)-»<l} 



(1-U) 



-20-4 



{(l-u)- e <£} 



V(l -u) 



{(l-u)-»<l} 



+ l) 2 (l-«) 



-20-4 



Vu < C. 



So we can let I — > oo and we get that (1 — u) 6 1 € 1 (SI) L 2 (f2), which means that f n ^ 1 _ u y-( 1+g) < C. 

By iterating the above argument for 6i + 1 = jfz^(0i-i + 1) for i > 1 and starting with 6*o = 0, we find 
that 1/(1 -tt) G L^Q) for all q G (l,oo). 

Standard regularity theory for elliptic problems applies again to give that 1/(1 — u) G C 2 ' Q (S1). Therefore, 
u is a classical solution and there exists a constant < K(C,N) < 1 such that || u || c(n) < K(C.N) < 1. 
This completes the proof of Theorem 3.6. ■ 



Theorem 3.7. For any dimension 1 < N < 8, f/iere exists a constant < C(N) < 1 independent of X such 
that for any < A < A*, the minimal solution u\ satisfies \\ u x || c(ri) < C(N). 

The theorem, which gives Theorem 1.2(2), will follow from the following uniform energy estimate on the 
minimal solutions u x . 

Lemma 3.8. There exists a constant C{p) > such that for each A G (0, A*), the minimal solution u x 
satis 



isfies || (i-l )d IUp(o) < C{p) as long as p < 1 + § + : 



w 2 dx < — I wAwdx = / |Vw| 2 dx, 



Proof: Proposition 3.3 implies that 

f 2f(x) 

L(i-%) 3 

for all < A < A* and nonnegative w G Hq(Q). Setting 

w = (1 - u A )' - 1 > , w/iere - 2 - V6 < i < , 

then H3.14f) becomes 

2[l-(l-«J*] a /(») 



: / (l- U J 2 - 2 |V W J 2 dx>A / 



(1-"J S 



dx . 



(3.14) 



(3.15) 



(3.16) 



On the other hand, multiplying (S)\ by jz^[(l — u x ) 21 1 — 1] an d applying integration by parts yield that 



i' / (1 — u A ) | Vu A \ dx — A 



2« - 1 



And hence l|3.16|l and (|3.17|l reduce to 

a* 2 r f(x) 



2i-U n (1-wJ 2 

,'2 



> A(2 



da; - 2A 



/(*) 



/(*) 



(1-«J 3 
-dx . 



dx + 4A 



/(*) 



(l-uJ»-< 



dx . 



dx 



(3.17) 



(3.18) 



2i-lVn (1-«J 3 - 24 
From the choice of i in (|3.15|l we have 2 + > 0. So l|3.18|l implies that 



L (i 



/(*) 



\3-2j 



< C 

< c 

< C 



3-i 



3-i 
/3-2« 



/(*) 



dx 



dx 



(3.19) 



dx 



where Holder's inequality is applied. From the above we deduce that 

/(*) 



(i-« x ) 3 -» 



dx < C. 



(3.20) 



Further wc have 



/(*) 



(1-«J 3 



2 i 

3 



f =™ f 

Jn (1-«J 3 (321) 

r / T, -r^pdx <C. 

(l-u x ) 3 - a « - 



Therefore, we get that 



where -in view of (|3.15|) - 



fix) 



<C, (3.22a) 



p=l_J!< 1 + 1 + 2,/H. (3.22b) 



Proof of Theorem 3.7: This follows from Lemma 3.8 and Theorem 3.6, where p — 1 when the dimension 

f<l + t+V§ 



N = 1, p can be taken to be 1 + 1 when N — 2. For N > 2, the reasoning applies as long as^ < 1 + 1 + 2^ ' 2 



which happens when N < 8. ■ 
Finally, we note the following easy comparison results and we omit the details. 

Corollary 3.9. Suppose fi,fs ■ ft — > (0,1] satisfy fi(x) < fz(x) on ft, </ien A*(ft,/i) > A*(f2,/2) and for 
< A < A*(f2, /i) we /icrae iti(A,a:) < M2(A,cc) on Q, where u\(X,x) (resp., U2{X,x)) are the unique minimal 
positive solution of 

— Au = ^/j;^ (resp., —Am = -p-^p-) on anrf u = on dfl. 

Moreover, if fi(x) > f2(x) on a subset of positive measure, then ui(A, x) < Ua(A, x) for all x £ Q. 

We note that if one considers the cases of power-law or exponential profiles for (S)\ defined in a ball, 
then the minimal positive solution corresponds to the lowest branch in the bifurcation diagram, the one 
connecting the origin point A = to the first fold at A = A*, see section §5. 

3.3 Existence of solutions at A = A* 

In this subsection, we study the existence of positive solutions at the critical voltage A = A*. We first deal 
with the existence of minimal solutions m a * for (S')a*- 

Lemma 3.10. Suppose there exists < C < 1 such that \\ u x || c(fi) < C for each A < A* . Then u x , = 
lim A/A „ u x exists in the C 2 ' Q (0) topology for some < a < 1. Moreover, there exists 6 > such that 
the solutions of (S)\ near (A*,u A ») form a curve p(s) = {(A(s), v(s)) : |s| < 5}, and the pair (X(s),v(s)) 
satisfies: 

A(0) = A*, A'(0) = 0, A"(0) < 0, and «(0) = u A „, v'(0)(x) > in ft. (3.23) 

Proof: The proof is similar to a related result of Crandall and Rabinowitz cf. [H] [7], so we will be brief. 
Firstly, the assumed upper bound on u x in C 1 and standard regularity theory, show that if / G C(ft) then 
II u x \\ c2 a(Ci) < C for some < a < 1 (while if / G L°°, then || u x \\ cl>a (fl < C). It follows that {(\,u x )} is 

precompact in the space K x C 2 ' a , and hence we have a limiting point (A*, u x , ) as desired. Since is 

nonnegative, Theorem 3.2 of [Hj characterizes the solution set of {S)\ near (\* ,u x , ): A(0) = A*, A'(0) = 0, 
v(0) = u x , and v'(Q) > in ft. Finally, the same computation as in Theorem 4.8 in |B] gives that A"(0) < 0. 
■ 

Remark 3.3. Lemma 3.10 implies that if the minimal solution u x (x) satisfies || u x || . fl . < C < 1 (which 
occurs when N < 8), then there exists two distinct solutions for (S)\ for A in a deleted left neighborhood of 
A* . A version of this result will be established variationally in an upcoming paper. 



The following theorem gives the uniqueness of (classic) solutions for (S)\* 



Theorem 3.11. Suppose there exists < C < 1 such that || u x || C(S) < C for each A < A*. Then the minimal 
solution u x , = lim A/ , A „ u x obtained above satisfies the following properties: 

1. The smallest eigenvalue ju(A) at A = A* of the linearized operator £« a ,a = —A — (^Tp 071 ^ * s zero - 



2. u x , is the unique solution of (S) 



A* • 



Proof: (1) Applying Proposition 3.3(2) we sec that /i(A) > on the minimal branch for any A < A*, hence 
the limit /i(A*) > 0. If now fJ.(X*) > the Implicit Function Theorem could be applied to the operator 
Lu x , ,a* , and would allow the continuation of the minimal branch A i— » u x of classical solutions beyond A* , 
which is a contradiction and hence = 0. 

(2) Suppose now u is any solution such that u > u x , . Since ^(A*) = 0, let </> be any positive eigenfunction 
in the kernel of L Ux , t \* and write, 

which yields that 

On the other hand, since — Ac/) — ^jz^f^p <f>, we have 



Since the integrand is nonnegative it follows that 



112 

+ - —(u-u xt ) a.e.inn. (3.24) 



{i-uf (i-u x ,r (i-«„) 



3 ' 

A* 



If now || u || i0 o< C < 1, then u is a classical solution as in Theorem 3.6, and we conclude that u = u x , on 
Q. U 

Now Theorem 1.3 is a direct result of Theorem 1.2(2) (or Theorem 3.7), Lemma 3.10 and Theorem 3.11. 



4 Uniqueness and Multiplicity of Solutions 

The purpose of this section is to discuss uniqueness and multiplicity of solutions for (S)\. Note that Lemma 
3.10 gives that for some < A2 < A*, there exists at least two solutions for (S)\ with A <G (A£, A*), which is 
Theorem 1.4(2). In the following we shall focus on the uniqueness when A is small enough. We first define 
non-minimal solutions for (S)\ as follows: 

Definition 4.1. A solution < u < 1 is said to be a non-minimal positive solution of (S)\, if there exists 
another positive solution v of (S)\ and a point x £ Q such that u(x) > v(x). 

Lemma 4.1. Suppose u is a non-minimal solution of (S)\ with A e (0, A*). Then the smallest eigenvalue 
yu(A) of the linearized operator L u ^\ = —A — on Q must be negative. 

Proof: For any fixed A G (0, A*), let u x be the minimal solution of (S)\. Clearly we have w = u — u x > 
in O, and 

-Aw- 1 \ 2 W = ° inQ - 

Hence we deduce from the strong maximum principle that u x < u in fi. 



Let n,Q — {x £ : f(x) = 0} and f2/f2o = {x £ fl : f(x) > 0}. Direct calculations give that 

2A/ 1 1 2 Jo, xe^o; 

(i-") 3 - /L (i- U ) 2 (i-« x ) a (i-«) 3V A ' J | <0 , xe n/n . 



(4.2) 



From this we get 

A ivn D /[ (T^F - (T^F - (T^F (U - (« - «J < o. 

Now suppose that /x(A) > 0. Then for each 4> £ Hq(Q) we have 

(L u , x <f>, </>)=[ (| V0| 2 - ^%d> 2 ) > . (4.3) 
Putting <f> = u — it x in (|4.3|) , we get from the left equality of jO} that 

A / n/O0 /[ (T^F " (T^F " JTW {u " Mj] (u " " J 
= A // [ 0^F " " (I^F (u " uJ] (u ~ uJ - 

which contradicts l|4.2(l . and we are done. ■ 

Remark 4.1. Proposition 3.3 and Lemma 4.1 give that for < A < A*, the smallest eigenvalue of L m ,a is 
necessarily negative if u is a non- minimal solution of (S)\, while the smallest eigenvalue of L u ^\ is positive 
if u = u x is a minimal solution of (S)\. For parabolic problems of the type (1.1), it is well-known that the 
spectrum of the linearized operator about any steady-state solution determines the stability of solutions for 
(1.1). Therefore, u x (x) is the unique stable steady-state of (1.1). In our upcoming paper [TT], we shall prove 
that the dynamic solution of (1.1) with A < A* (and A = A* for N < 8) will globally converge to its unique 
minimal solution u x (x). 

Now we are able to prove the following uniqueness result, which completes the proof of Theorem 1.4. 

Theorem 4.2. For every M > there exists < \\{M) < A* such that for A £ (0,X*(M)) the equation 
(S)\ has a unique solution v satisfying: 

1- II hi \3 111 ^ M as long as the dimension N = 1. 



(1 


-v) a 




f 


(1 




/ 1 


1 

„\3 



3- 11(1^11^/2 <M andN>2. 

Proof: For any fixed X £ (0, A*), let u x be the minimal solution of {S)\ and suppose (S)\ has a non-minimal 
solution u. Lemma 4.2 then gives 

,_. >,o , f 2A(w - u x ) 2 f(x) , 

\V(u-u.)\ 2 dx < / — y — XJ J y ' dx. 

Jn {1-uf 

This implies in the case where N > 2 that 



C{N)(J^{u-u x )^dx) * <xj^^L{u-u x ?dx 



< 2A 



/ 



(1-u) 3 
<2XM&( (u- u x )-^dx^j 



) N (J^(u-u x )"-idx) 



N -2 
N 



which is a contradiction if A < — unless u = u.. If N = 1. then we write 

2MW 



C(1)||(u-«J||2 <a/ -^%(7i- UA ) 2 dx<2A||( w -u A )||L / Tr-^dx 
Jo, I 1 — M J Jo l 1 — U J 

and the proof follows. A similar proof holds for dimension N = 2. ■ 

Remark 4.2. The above gives uniqueness for small A among all solutions that cither stay away from 1 or 
those that approach it slowly. We do not know whether if A is small enough, any positive solution v of (S)\ 
satisfy / n (l — v)^^dx < M for some uniform bound M independent of A. Numerical computations do show 
that we may have uniqueness for small A -at least for radially symmetric solutions- as long as N > 2. 



5 Steady- State: Case of Power-Law Profile 

The issues of uniqueness and multiplicity of solutions for (S)\ with < A < A*, and even mere existence 
for (S)\* with N > 8 seem to be quite challenging problems. In this section, we discuss these problems in 
the case where / has a power-law permittivity profile, i.e, f(x) — \x\ a (a > 0). We shall also consider the 
domain to be a unit ball -Bi(O) C M w (N > 1) and A S (0, A*]. In this special case, the solutions of (S)\ 
must be radially symmetric, and (S)\ is then reduced to the following problem 

N — 1 \r a 

-U rr U r = —, 0<T<1, 

r {1-uY (5.1) 

u'(0)=0, u(l)=0. 

Here r = \x\ and < u = u(r) < 1 for < r < 1. 
Looking first for a solution of the form 

u(r) = 1 — (3w(P) with P = jr , 

where 7, (3 > 0, equation (|5.1|) implies that 

2 01 n , N - 1 XP a 1 

7 p(w H — w ) = 



P ' (3 2 j a w 2 ' 
We set w(0) = 1 and A = j 2+a /3 3 . This yields the following initial value problem 

., N - 1 , P a 



w 



-w' = — =- , P > . 



P w 2 (5.2) 

io'(0) = 0, io(0) = l. 



Since = we have (3 — l/w(j). Therefore, we conclude that 

1 



u(0) = 1 - 
7 2 



(5.3) 

A 



W 3 ("f) ' 

where ^(7) is a solution of H5.2J) . 

As done in one can numerically integrate the initial value problem (|5.2|l and use the results to 
compute the complete bifurcation diagram for (|5.1J) . We show such a computation of u(0) versus A defined 
in H5.3fl for the slab domain (N = 1) in Fig. 2. In this case, one observes from the numerical results that 
when N = 1, 

• There exists a unique solution for (S)\* ; 

• For < a < 1, there exist exactly two solutions for (S')a whenever A € (0, A*); 



Figure 2: Plots of \u(0)\ versus A for the power-law permittivity profile f(x) = \x\ a (a > 0) defined in the 
slab domain (N — 1). The numerical experiments show that it seems to exist a constant a* > 1 (analytically 
given in (|5.#H ). such that the bifurcation diagrams are greatly different for different ranges of a: < a < 1, 
1 < a < a* and a > a* . 



• For a > 1, it is however difficult to see in any other case the bifurcation diagram as u(0) — > 1. 

This leads us to the question of determining the asymptotic behavior of w(P) as P — > oo. Towards this 
end, we proceed it as follows. 

Set T) = logP, w(P) — P B V(rj) > for some positive constant B. Then we obtain from (|5.2f> that 

pa-2B 

pB-iyii + (2B + N - 2)P B ~ 2 V' + B(B + N - 2)P B ~ 2 V = y2 . (5.4) 
Choosing B — 2 = a — 2B so that B = (2 + a)/3, we get that 

V" + 3iV+ 3 2Q -V + (2 + Q)(3 ^ + Q - 4 V = i, . (5.5) 

We notice from (|5.5|) that only for the case where N > 2 or N — 1 with a > 1, that the equilibrium point 
V e of (|5.5|l must be positive and satisfies 

^= (2 + a )(3^ + a -4) >0 - ^ 

When N = 1, this is consistent with the numerical observation of Fig. 2. Linearizing around this equilibrium 
point by writing 

V = V e + Ce ar ' 0<C«1, 



we obtain that 



This reduces to 



2 3iV + 2a-2 (2 + a)(37V + a-4) 
, + § * + =0. 

3JV + 2a-2 , VS 
ff± = g ± — , (5.7a) 



with 



A = -8a 2 - (247V - 16)q + (97V 2 - 84iV + 100) . (5.7b) 
We note that a± < whenever A > 0. Define now 



Next, we discuss on N and a by considering the sign of A: 
Case (1). N and a satisfy either one of the followings: 

N = 1 with 1 < a < a* 



N > 8 with < a < a** 



(5.9a) 
(5.9b) 



In this case, we have A > and 

9 



V 



.(2 + a)(3iV + a-4) 

Further, we conclude that 

D 2+« / 9 
w ~ P 3 



Cie e ^ 



• • • , as 77 — > +cxd . 



,(2 + a)(37V + a-4) 



N-2 1 y7? 

CiP — +— h , as P^+oo, 



Since A = 7~ +Q /w (7), we have 



. , (2 + a)(3iV + a-4) 

A ~ A* = as 7 — > 00 . 



(5.10) 



For this case, we compute the numerical results for H5.9a(l in Fig. 2 and for (|5.9b|l in Fig. 3(b), respectively. 
Case (2). N and a satisfy any one of the following three: 



N = 1 with a > a* ; 

2 < N < 7 with a > : 
N > 8 unift. a > a** . 



(5.11a) 

(5.11b) 
(5.11c) 



Number of solutions for the case f(x) = \x 





A = Ai(< A*) 


A < A* 


A = A* 


< a < 1 




2 


1 


1 < a < a* 


1 


> 1 


1 


a > a* 


00 


> 1 


1 



and N = 1 



Table 3: Number of solutions to (1.1) which is defined in a unit ball Pi(0) with N = 1, where Ai 
(2+a)(a-i) ^ a * _ _i _|_ ^^27/2 and /(a;) = |a;| a is chosen to be power-law permittivity profile. 



Number of solutions for the case f(x) = |x| Q and 2 < N < 7 





A = A*(< A*) 


A < A* 


A = A* 


a>0 


00 


> 1 


1 



Table 4: Number of solutions to (1.1) which is defined in a unit ball Pi(0) with 2 < N < 7, where A* 
(2+a)(3A r +a 4) ^ f (x) = \x\ a is chosen to be power-law permittivity profile. 




x. x' x x' x . x' X 



Figure 3: Left figure: Plots of \u(0)\ versus A for the power-law permittivity profile f(x) = \x\ a (a > 0) 
defined in the unit ball P>i(0) C R N with 2 < N < 7. In this case, |u(0)| oscillates around the value A* 
defined in $5.10} and there exists a unique solution for (S)\* . Right figure: Plots of |u(0)| versus A for the 
power-law permittivity profile f(x) — \x\ a (a > 0) defined in the unit ball -Bi(O) C with N > 8. The 
characters of the bifurcation diagrams depend on different ranges of a: when < a < a**, there exists a 
unique solution for (S)\ with A € (0, A*) and there does not exist any solution for (S)\*; when a > a**, 
|u(0)| oscillates around the value A* defined in Q\5.1U\i and there exists a unique solution for (S)\*. 




Figure 4: Plots of |u(0)| versus A for the constant permittivity profile f(x) = 1 defined in the unit ball 
Bi(0) C with different ranges of N . In the case ofN>8, we have A* = (6 AT - 8)/9. 



Number of solutions for the case f(x) = \x\ a and N > 8 





A = A*(< A*) 


A < A* 


A = A*(= A*) 


A = A* 


< a < a** 




1 








a > a** 


oo 


> 1 




1 



Table 5: Number of solutions to (1.1) which is defined in a unit ball -Bi(O) with N > 8, where A* = 
(2+a)(3N+a-4) ^_ ^ q < q- < a** ), a** — 4 ~ 6N + 3 ^( N ~ 2 ) an( j. f (x) = \x\ a is chosen to be power-law 
permittivity profile. 

Number of solutions for the case f(x) = 1 





A = A*(< A*) 


A < A* 


A = A*(= A*) 


A = A* 


N = 1 




2 




1 


2 < N < 7 


oo 


> 1 




1 


N>8 




1 









Table 6: Number of solutions to (1.1) which is defined in a unit ball B>i(0), where A* = 6N 9 8 and f(x) = 1 
is chosen to be constant permittivity profile. We note that A* = A* = (6N — 8)/9 for N > 8. 



In this case, we have A < and 
9 



V 



{2 + a){3N + 
Further, we conclude that 
9 



4) 



w 



2 + a 
P — 



■Cie 



■CiP- 



-A 



'cos 



N-2 / \ j — A 

— cos(- logP + C" 2 



.(2 + a)(3N + a-4) 

And we also obtain from A = ^ 2+a / 'w 3 (7) that 

. _ (2 + a)(3N + a - 4) 
A ^ A* — 



as 77 



P 



+00 . 



-00 . 



(5.12) 



7 



When N and a separately satisfy (|5.11a(l . (|5.11b(l and l|5.11c(l . the typical diagrams are computed in Fig. 2, 
Fig. 3(a) and Fig. 3(b), respectively. 

Combining the numerical results of Fig. 2 and Fig. 3, we plot the bifurcation diagram of the constant 
permittivity profile defined in the unit ball with different ranges of N . The result of such a computation is 
shown in Fig. 4, from which one can observe the uniqueness, (infinite) multiplicity of the solutions for (S)\ 
with A e (0, A*) and different ranges of N: N = 1, 2 < N < 7 and N > 8. We note that A* = (6 AT - 8)/9 
when N > 8. 

Applying above numerical results, in Table 3~5 we give the number of solutions for H5.1J) depending on 
and a. In Table 6, we give the number of solutions for (|5.1|l with constant profile f(x) = 1, which shows 
that N = 8 is the critical dimension for (S3) - 

Remark 5.1. Under the assumptions of (|5.11J1 . since ui(P) in H5.12(l is oscillatory for P >> 1, we expect that 
|u(0)| oscillates around the value A* = ( 2 +")( 3 ^ , + a ~" 4 ) as p _> j n particular, this implies that for this 
case, (S)\ has infinitely multiple solutions for (S)\ t . We compute the numerical results of this case in Fig. 2 
and Fig. 3, from which we can observe the uniqueness, (infinite) multiplicity of the solutions for (S)y. when 
N and a satisfy any one of the three cases in H5.11fl . then there exists a series of {A^} satisfying 



A = 0, 
Ai = A* 



A 2 fe / K a 
A2/C-1 \ A* 



00 : 



as 



such that: there exist exactly 2k + 1 solutions for (S)\ with A S (A2fe, A2&+2); and there exist exactly 
2k solutions for (S)\ with A G (\2k+i , ^2k-i)', further, there exist infinitely multiple solutions for (S')a». 



Therefore, it is reasonable to believe that the multiplicity of solutions for the general (S) \ greatly depends 
on the permittivity profile f(x), the dimension TV and the value of A. 

Remark 5.2. Our results show that for f(x) = \x\ a with N > 8 and < a < a** , then there does not exist 
any classic solution for (S)\*, where A* = (2 + a)(3N + a — 4)/9; but for other cases of N and a, there 
exists a unique solution for (S)\* . Therefore, for N > 8 it seems from these results that whether there exist 
solutions for (S)\* depends on the varying permittivity profile f(x). However, we conjecture that for N > 8 
there is no solution for (S) \* if the permittivity profile f(x) > on Q. 
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